1/1 


AD-A161  549  MULTIPLE-USER  INTERFERENCE  IN  A  FREQUENCY -HOPPED 

SPREAD-SPECTRUM  COMMUNIC  (U)  ILLINOIS  UNIV  AT  UR8ANA 
COMPUTER  SYSTEMS  GROUP  A  U  LAM  SEP  84  R-1816 
UNCLASSIFIED  N88014-84-C-0149  F/G  17/2  NL 


4$  ’% 


REPORT  R- 1016  SEPTEMBER,  1984 


Ul  LU-ENG  2210 


&  s COORDINATED  SCIENCE  LABORATORY 


MULTIPLE 'USER  INTERFERENCE 
IN  A  FREQUENCY-HOPPED 
SPREAD -SPECTRUM  COMMUNICATION 
SYSTEM  USING  DIFFERENTIAL 
PHASE -SHIFT  KEYING 


UNIVERSITY  OF  ILLINOIS  AT  URBANA-CHAMPAIGN 

**  "  18  8» 


Unclassified 


security  classification  of  this  page 


la.  report  SECURITY  classification 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OECLASSIFICATION/OOWNGRAOING  SCHEDULE 

N/A  _ 


4.  PERFORMING  ORGANIZATION  REPORT  NUMSERIS) 

R-1C16  (TOLD- ENG  84-2210) 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 

None 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release,  distibution 
unlimited 


8.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

N/A 


fla.  NAME  OF  PERFORMING  ORGANIZATION  Bh.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

Coordinated  Science  Laboratory  j  dfwueabitt 


University  of  Illinois 


6c.  AOORESS  (City.  Stata  and  ZIP  Coda) 

1101  W.  Springfield 
Urbana,  IL  61801 


Ba.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 


Joint  Services  Electronics  Program 


7b.  AOORESS  (City,  Stau  and  ZIP  Coda) 

Arlington,  VA  22217 


Bb.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(if  applicable) 


N00014-84-C-0149 


Be.  AOORESS  (City,  State  and  ZIP  Code) 

10.  SOURCE  OF  FUNDING  NOS. 

Arlington,  VA  22217 

PROGRAM 
ELEMENT  NO. 

PROJECT 

NO. 

TASK 

NO. 

WORK  UNIT 
NO. 

11.  TITLE  (include  Security  Classification) 

Multiple-User  Interference  in  a  ... 

N/A 

N/A 

N/A 

N/A 

12.  personal  authorisi 
Lam,  Alex  W. 


13a.  type  of  report 

Technical 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVEREO 
FROM _  TO 


14.  DATE  OF  REPORT  (Yr..  Mo..  Day)  15.  PAGE  COUNT 

August  1984  63 


17. 


FIELD 


COSATI  COOES 


18.  SUBJECT  TERMS  (Continue  on  reverie  if  neceseary  and  identify  by  block  number) 

Near— '^timum  FH-DPSK  receiver,  Reed-Solomon  frequency  hop¬ 
ping  patterns ..multiple-user  interference .error  probabilities 

\  \ 


19.  ABSTRACT  (Continue  on  revene  if  neceaeary  and  identify  by  block  number (' 

A  design  of  short  Reed-Solomon  hopping  sequences  which  have  excellent  correlation 
properties  is  proposed.  Within  the  design,  any  two  sequences  will  have  at  most  one 
partial  hit  or  one  full  hit  when  used  asynchronously.  The  behavior  of  a  near-optimum 
differential  phase-shift  keying(DPSK)  receiver  with  multiple-user  interference  over  a  non¬ 
fading  channel  is  discussed.  Bit  error  rates  for  the  receiver  with  single  interference 
in  one  chip  and  the  receiver  with  signle  interference  in  J  distinct  chips  are  found. 
Finally,  generalization  of  the  multiple-user  interference  model  is  considered. 


20.  OISTRIBUTION/AVAILABILITY  of  abstract 
unclassifieo/unlimited  2  same  as  rpt.  Z  otic  users  □ 


22*.  NAME  OF  RESPONSIBLE  INOIVIOUAL 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassif ied 


22b.  TELEPHONE  NUMBER 
(include  Area  Code) 


22c.  OFFICE  SYMBOL 

N/A 


ied 


MULTIPLE-USER  INTERFERENCE  IN  A  FREQUENCY-HOPPED  SPREAD-SPECTRUM 
COMMUNICATION  SYSTEM  USING  DIFFERENTIAL  PHASE-SHIFT  KEYING 


BY 

ALEX  WAIHO  LAM 

B.S.,  University  of  Illinois,  1982 


THESIS 

Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Electrical  Engineering 
in  the  Graduate  College  of  the 
University  of  Illinois  at  Urb ana -Champaign,  1984 


Urbana,  Illinois 


ft. 


Accession 

NTIS  0"  t\  i 
DTIC  T 
IKv 


Fv 


MULTIPLE-USER  INTERFERENCE  IN  A  FREQUENCY-HOPPED  SPREAD-SPECTRUM 
COMMUNICATION  SYSTEM  USING  DIFFERENTIAL  PHASE-SHIFT  KEYING 

by 

Alex  Waiho  Lam 

Department  of  Electrical  Engineering  and 
The  Coordinated  Science  Laboratory 
University  of  Illinois  at  Urbana-Champaign,  1984 

ABSTRACT 

A  design  of  short  Reed-Solomon  hopping  sequences  which  have  excellent 
correlation  properties  is  proposed.  Within  the  design,  any  two  sequences 
will  have  at  most  one  partial  hit  or  one  full  hit  when  used  asynchronously. 

The  behavior  of  a  near-optimum  differential  phase-shift  keying  (DPSK)  receiver 
with  multiple-user  interference  over  a  nonfading  channel  is  discussed.  Bit 
error  rates  for  the  receiver  with  single  interference  in  one  chip  and  the 
receiver  with  single  interference  in  J  distinct  chips  are  found.  Finally, 
generalization  of  the  multiple-user  interference  model  is  considered. 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  his  sincere  gratitude  to  his  advisor. 
Professor  Dilip  V.  Sarwate,  for  his  guidance  and  expert  assistance  during 
the  preparation  of  this  thesis.  The  author  also  wishes  to  thank  Mrs. 
Phyllis  A.  Young  for  her  excellent  typing  of  the  manuscript. 


V 


TABLE  OF  CONTENTS 


CHAPTER  Page 

I  INTRODUCTION .  1 

II  HOPPING  PATTERNS  FOR  THE  FREQUENCY-HOPPED 

MULTIPLE  ACCESS  COMMUNICATION  SYSTEM .  13 

2.1  A  Hopping  Pattern  Design  for  the 

Asynchronous  Upstream  Communications .  14 

2.2  A  Hopping  Pattern  Design  for  the  Synchronous 

Downstream  Communications .  23 


III  PROBABILITY  OF  ERROR  ANALYSIS  OF  A  NEAR-OPTIMUM  FH-DPSK 

RECEIVER  WITH  ONE  INTERFERENCE  AND  A  NONFADING  CHANNEL .  26 

IV  PROBABILITY  OF  ERROR  ANALYSIS  OF  A  NEAR-OPTIMUM  FH-DPSK 
RECEIVER  WITH  MULTIPLE-USER  INTERFERENCE  AND  A 


NONFADING  CHANNEL .  42 

V  CONCLUSION .  53 

APPENDIX  A  THE  ERROR  PROBABILITY  P  (u,x) .  54 

e 

REFERENCES .  6  3 


CHAPTER  I 


INTRODUCTION 

In  1977,  Cooper  and  Nettleton  proposed  a  frequency-hopping  multiple-access 
(FHMA)  system  for  cellular  mobile  radio  communications  [1],  [2],  The  details 
and  advantages  of  this  system  have  been  discussed  extensively  [3],  [4 ]  by 
previous  authors  and  shall  not  be  repeated  here.  Basically,  half  of  the 
allotted  mobile  radio  spectrum  is  used  for  base-to-mobile  transmission  and 
the  other  half  for  mobile-to-base .  The  data  stream  is  partitioned  into  groups 
of  n  bits,  and  each  group  is  encoded  into  2n  =  N  channel  bits  using  an  (N,n) 
orthogonal  code.  The  codewords  are  taken  to  be  the  rows  of  an  N  x  N  Hadamard 
matrix  [5]  with  elements  ±  1.  Each  of  the  N  bits  is  transmitted  using  a 
different  carrier  frequency,  according  to  the  user's  unique  frequency-hopping 
pattern,  or  address.  A  differentially  coherent  phase-shift  keying  (DPSK) 
technique  is  used  to  modulate  each  bit  onto  its  particular  carrier  frequency. 
The  basic  time-frequency  encoded  waveform  can  be  written  as  [3] 

_ N-l 

s(t)  =  i/Z?  E  pT  (t-kT  )cos(27T(fQ+a^Af H+G^)  (1.1) 

k=0  c  C 


where  P  denotes  signal  power 


PT  (t) 
c 


is  the  rectangular  pulse  function: 


0  <  t  <  T 
—  c 

otherwise 


(1.2) 


T  is  the  time-chip  duration,  9  denotes  the  phase  of  the  k-th  radio- 

C  K. 

frequency  (RF)  pulse  and  contains  the  binary  DPSK  information,  f  is  the 
nominal  carrier  frequency,  Af  is  the  minimum  frequency  shift,  and  a  is 


the  element  of  the  address  that  specifies  the  k-th  carrier  frequency  of  the 


hopping  pattern.  If  we  let  the  waveform  duration  be  T,  we  have  that  T  =  NT^ , 
where  N  is  the  number  of  chips  per  waveform.  Each  time-chip's  phase  is 
compared  to  the  corresponding  time-chip's  phase  in  the  previous  waveform  to 
obtain  the  DPSK  information. 

The  model  for  the  DPSK  receiver  used  in  [1]  and  [2]  is  shown  in  Fig.  1.1. 
Cooper  and  Nettleton  studied  the  performance  of  this  receiver  via  Monte 
Carlo  computer  simulation.  Martersteck  [3]  has  recently  derived  an  optimum 
receiver  for  the  signaling  format  just  discussed  (Eq.  (1.1))  and  has  also 
studied  a  simplified  version  of  the  receiver  which  is  nearly  optimal  for  small 
signal-to-noise  ratio  (SNR).  He  has  shown  that  the  bit  error  rate 
performance  of  the  near-optimum  receiver  is  better  than  that  of  the  Cooper- 
Nettleton  receiver,  in  both  nonfading  and  fading  environments  with  no 
interference.  An  outline  of  the  aforementioned  derivation  is  briefly  given 
below. 

A  receiver  for  the  DPSK  Cooper-Nettleton  system  (assuming  a  nonfading 
channel)  must  decide  on  the  transmitted  codeword  based  on  the  received 
waveform 

_ N-l 

r(t)  =  /iP  E  pT  (t-kTc)cos(wkt  +  0^ 
k=0  c 

N-l 

+  /2P  E  PT  (t-T-kTc)hikcos(mkt  +  0k)  +  n(t)  (1.3) 

k=0  c 


where  w  =  2ir(f  +a  Af),  h  is  the  ±1  element  in  the  i-th  row  and  k-th  column 

tC  w  1C  JlJx 


of  the  Hadamard  matrix,  n(t)  is  the  additive  white  Gaussian  noise  (AWGN)  due 
to  receiver  thermal  noise,  0  is  a  random  phase  variable,  and  (the  i-th  row 


4 


of  the  Hadamard  matrix)  is  the  transmitted  codeword.  Let  denote  the 
hypothesis  that  the  i-th  row  (s^)  of  the  Hadamard  matrix  is  being  transmitted. 


Defining 


(k+l)T 

X,  =  /  cr(t)cos  a)  t  dt, 

k  kT  k 

c 


(k+l)T 

Y  =  /  cr(t)sinm  t  dt 

R  kT  k 

c 


T+(k+l)T 

X,'  =  /  cr(t)cosu),  t  dt, 

*  T+kT  k 

c 


T+(k+l)T 

Y '  =  /  cr(t)s into,  t  d t  (1.4) 

k  T+kT  k 

c 


and  using  optimum  receiver  theory,  the  optimum  (nonfading)  receiver  will  base 
its  decision  on  the  statistic 


A(H.) 


N-l 

A  e"uo<(8PNo2[(\'H'ikx()2  +  'WP2”  >>•  i-°-1 


.N-l 

(1.5) 


where  Iq(-)  is  the  zero  order  modified  Bessel  function  of  the  first  kind.  The 

optimum  receiver  must  then  calculate  A(H^)  for  i  =  0,1,..., N-l  and  decide 

was  transmitted  whenever  A(H^)  >  A(H_.)  for  all  j  =  0,1,..., N-l  and  j  ^  i. 

In  other  words,  the  optimum  receiver  employs  a  maximum  likelihood  decision  rule. 

The  optimum  receiver  structure  can  be  simplified  by  approximating  the 

2 

function  JLnlg(x)  in  Eq.  (1.5).  In  particular,  inl^x)  x  /4  for  x  <<  1,  and 
tolg(x)  «=>  x  for  x  >>  1.  Thus,  for  small  bit  SNRs,  Eq.  (1.5)  becomes 


*  N-i  2 

A  (Hi)  =  Z  +  (Yk+h.kY^)  ]  ,  i  =  0,1,..., N-l 

k=0 


(1.6) 


which  further  reduces  to 


* 

A  (H.) 


N-l 


k=0 


<V£  +  \Yk)hik 


,  i  =  0,1 , . . . ,N-1 


(1.7) 


5 


Either  of  these  two  expressions  determines  the  structure  of  the  so-called 

k 

small-SNR  receiver.  The  former  expression  for  A  (H_^)  is  easiest  to  work  with 
analytically,  whereas  the  latter  expression  is  easiest  to  implement  since 
fast  Hadamard  transforms  may  be  used  [4].  Curiously  enough,  the  small-SNR 
receiver  can  also  be  shown  to  be  optimum  in  a  Rayleigh-fading  environment 
provided  that  and  Y^,  k  =  0,1,..., N-l  are  uncorrelated.  For  the  scheme 
considered  by  Martersteck,  the  X  and  Y  are  nearly  uncorrelated.  Thus,  the 

K.  K 

small-SNR  receiver  is  nearly  optimum  over  a  fading  channel.  We  shall  adopt 
the  convention  of  [ 3]  and  refer  to  the  small-SNR  receiver  as  the  near -optimum 
receiver.  A  block  diagram  of  this  receiver  is  shown  in  Fig.  1.2. 

For  large  chip  SNRs,  Eq.  (1.5)  yields 


**  a 

A  (H.)  =  Z  [(Xk+hikXkr  +  (Wk}  1  ’  1  =  °’1 . N_1  (1*8) 

k=0  '  l" 


and  we  call  the  receiver  induced  by  this  expression  the  large-SNR  receiver. 

The  upper  bound  (union  bound)  for  the  bit  error  rate  performance  was 
obtained  analytically  by  Martersteck  for  the  near-optimum  receiver  over  a 
fading  and  nonfading  channel  with  no  interference.  Recently,  McClatchey  [4] 
has  found  that  the  bit  error  probability  for  the  large-SNR  receiver  and  the 
near-optimum  receiver  (for  a  nonfading  channel)  are  both  upper  bounded  by  the 
same  function.  Since  the  union  bound  is  known  to  be  tight  for  large  SNRs, 
the  two  receivers  thus  have  nearly  the  same  performance  at  large  SNRs. 
Furthermore,  we  would  expect  the  near-optimum  receiver  performance  to  be  better 


.  - 
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than  the  large-SNR  receiver  performance  (for  a  nonfading  channel)  at  small 
SNRs,  since  the  former  receiver  is  based  on  a  small— SNR  approximation  whereas 
the  latter  is  based  on  a  large— SNR  approximation.  Finally,  we  know  that  the 
near-optimum  receiver  is  nearly  optimum  in  a  Rayleigh-f ading  channel.  It 
follows  that  no  further  consideration  of  the  large-SNR  receiver  is  necessary. 

Figure  1.3  (from  [4])  shows  the  simulation  data  (via  Monte  Carlo  simula¬ 
tions)  for  the  optimum  and  near-optimum  receiver  for  N  =  32.  Upper-bound 
performance  for  the  near-optimum  receiver  and  simulation  performance  for  the 
Cooper-Net tleton  receiver  are  also  shown  in  Fig.  1.3.  Note  that  the  near¬ 
optimum  receiver  performance  (simulated)  is  better  than  the  Cooper-Nettleton 
receiver  performance  for  SNRs  above  3  dB  in  a  nonfading  no-interference 

environment.  In  a  digital  voice  communication  system,  a  bit  error  rate  of 
-3 

about  10  is  adequate  for  maintaining  reliable  communications.  At  N  =  32, 
this  corresponds  to  a  bit  SNR  of  about  8  to  10  dB.  Since  the  union  bound  is 
tight  for  large  bit  SNRs,  the  upper-bound  performance  of  the  near-optimum 
receiver  for  normal  usable  ranges  of  bit  SNR  is  thus  expected  to  represent 
the  receiver's  true  performance.  This  fact  may  be  verified  by  Fig.  1.3  for 
N  =  32  with  a  nonfading  no-interference  environment.  In  this  thesis,  it  is 
assumed  that  the  upper-bound  bit  error  rate  is  a  true  performance  measure  for 
the  near-optimum  receiver  when  a  sufficiently  large  bit  SNR  is  used. 

From  all  the  facts  that  we  have  seen,  the  near-optimum  receiver  is  by 
far  the  best  receiver  known  for  a  Cooper-Nettleton  system:  it  is  relatively 
simple  to  implement,  easy  to  synthesize,  is  analytically  tractable,  and  has 
performance  that  is  superior  to  the  Cooper-Nettleton  receiver  for  large  bit 
SNRs  in  a  no-interference  environment. 


J 

t  J* 


This  thesis  investigates  one  important  system  aspect  which  has  not  yet 

been  considered:  the  multiple-user  interference  problem.  In  a  FHMA  system, 

users  usually  share  the  same  frequency  band,  but  their  hopping  sequences  are 

not  necessarily  disjoint.  Consequently,  cross-user  interference  is  bound 

to  arise  if  large  numbers  of  mobiles  are  transmitted  simultaneously  in  the 

same  cell.  This  multiple-user  interference  can  degrade  FHMA  system  performance 

rather  seriously.  Indeed,  Yue  [6]  has  estimated  that  the  mobile-to-base 

number  of  users  per  cell  for  a  DPSK/FHMA  system  with  Cooper-Nettleton  receivers 

can  be  as  low  as  26,  assuming  random  addressing  and  bit  error  rates  less  than 
-3 

10  .  The  number  26  is  so  low  maybe  because  users  are  hopping  randomly  on  the 

spectrum  band.  In  this  thesis,  instead  of  random  hopping,  we  assume  that 
each  user  has  its  own  specific  nonrandom  frequency-hopping  sequence. 

If  the  mobiles  can  transmit  in  a  fashion  such  that  the  probability  of 
having  two  or  more  mobiles  transmitting  the  same  frequency  at  the  same  time  is 
small,  then  we  may  be  able  to  keep  the  multiple-user  interference  below  a 
tolerable  level.  In  the  analysis  of  [1},  the  multiple-user  interference  was 
modeled  as  a  random  white  Gaussian  process.  We  do  not  make  the  same  assump¬ 
tion  in  this  thesis.  However,  we  assume  that  signals  from  all  the  mobile 
transmitters  received  by  the  receiver  do  have  the  same  power. 

From  the  foregoing  content  of  this  chapter,  we  know  that  each  mobile  is 
transmitting  ±  1  channel  bits  (elements  of  the  rows  of  a  N  x  N  Hadamard 
matrix)  using  frequencies  according  to  its  own  unique  frequency-hopping 
pattern.  For  N  =  4,  the  Hadamard  matrix  is  simply 


>;s'v .  :mm-  <  vv 


r 


-4 


L 


10 


H  = 


(1.9) 


1111 
1-1  1-1 
1  1-1-1 
1-1-1  1 

Suppose  there  are  only  two  mobiles,  say  A  and  B,  transmitting  in  the  system. 


Without  loss  of  generality,  we  may  assume  that  A  is  transmitting  s^,  the 


0-th  row  of  H  (all  ones),  in  the  time  slot  [0,T].  If  a  receiver  A  at  the  base 


is  decoding  the  message  sent  by  A,  A  has  to  tune  at  A’s  transmitting 


frequencies  from  time  to  time  using  A's  unique  hopping  sequence.  The  receiver 


ft  ft  ft 

A  then  computes  A  (H^)  and  chooses  the  hypothesis  for  which  A  (H.)  is 


largest.  However,  in  order  to  compare  A  (H^)  and  A  (H^),  it  need  only  compare 


(X1X1  +  Y1YP  +  (X3X3  +  Y3Y3}  With  '(XIX[  +  YlYl}  “  (X3X3  +  Y3Y3}  becauSe 


h00  *  h10  and  \)2  =  h12*  N°W  consic*er  tbe  signal  of  B.  Suppose  that  B's 


signal  is  such  that  A  and  B  are  using  the  same  frequency  in  the  2nd  time-chip. 


However,  note  that  in  deciding  whether  A  (H^)  or  A  (H^)  is  larger,  the 


quantities  Y2^2  are  not  even  considered.  Thus,  the  presence  of  the 

interferer  makes  no  difference  in  this  decision.  As  a  second  case,  consider 


if  the  interference  from  B  occurred  in  the  1st  (or  3rd)  time-chip,  then  the 


*  * 

comparison  between  A  (H1 )  and  A  (H  )  will  inevitably  be  affected;  but  B  may 


V  “  “  '0' 

be  transmitting  a  +1  or  a  -1  bit.  If  B  is  sending  a  +1  bit,  the  interference 


will  be  constructive.  In  other  words,  interference  from  B  is  actually  helping 
* 


A  to  decide  correctly  that  A  is  transmitting  a  +1  bit,  so  the  error  probabil¬ 
ity  decreases.  On  the  other  hand,  if  B  is  sending  a  -1  bit,  the  interference 


will  be  destructive,  and  the  error  probability  increases.  As  a  summary,  in 


*  * 

computing  the  probability  P(A  (H. )  >  A  (Hn)),  an  interference  at  the  0-th 


v  0 

or  2nd  time  chips  will  have  no  effect  on  the  probability,  while  an 


^  -3  ^  .-3  ^  .O  *-"■  .W. 


>y- 
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interference  at  the  1st  or  3rd  time-chips  will  be  equally  likely  to  be 

constructive  or  destructive.  It  is  interesting  to  note  that  the  system  has 

a  sort  of  'built-in'  anti-interference  capability.  This  capability  is  a 

direct  consequence  of  the  system's  orthogonal  channel  coding  scheme,  and  is 

expected  to  increase  with  the  codeword  length  N. 

*  * 

Likewise,  we  can  compute  P(A  (H.)  >  A  (H  ))  for  all  i  and  j.  It  turns 

J  i 

out  that  this  probability  does  not  depend  on  j  or  i.  This  is  intuitively 
satisfactory  because  every  codeword  is  equally  likely  to  be  transmitted,  and 
the  interference  signal  is  randomly  distributed  among  the  time -chips.  The 

word  error  probability  can  then  be  upper  bounded  by  the  union  bound,  which 

A  k 

is  simply  3-P(A  (H^)  >  A  (H^))  in  this  case. 

For  multiple  interference  at  multiple  time-chip  periods,  the  analysis  is 

« 

more  involved,  but  it  still  follows  a  similar  line  of  reasoning.  In  most 
cases,  part  of  the  interference  signals  will  be  rejected.  Among  the  rest  of 
the  interferences  which  are  considered  by  the  receiver ,  some  of  them  will  be 


constructive,  while  the  others  will  be  destructive.  We  have  shown  that  the 

b  k  k 

word  error  probability  is  bounded  by  the  union  bound:  (N-l)P(A  (H.)  >  A  (H.)). 

.1  i 

We  will  also  show  that  the  bit  error  probability  is  bounded  by 


7  P(A*  (H  )  >  /(H)). 
-  J  l 


The  detailed  development  of  this  function  is  largely 


algebraic  and  manipulative  in  nature,  and  is  therefore  relegated  to  the 


Appendix. 


The  remainder  of  this  thesis  is  organized  as  follows.  In  Chapter  II  we 


discuss  a  design  of  good  frequency  hopping  sequences  in  which  any  two  sequences 
will  have  very  little  interference  over  all  time  shifts  asynchronously.  In 
Chapter  III  we  study  the  performance  of  the  near-optimum  receiver  with  only 


one  interfering  signal  at  a  particular  time-chip  period.  In  Chapter  IV  we 

present  an  upper  bound  of  the  bit  error  rate  for  the  near-optimum  receiver 

with  multiple-user  interference  assuming  that  each  time-chip  has  at  most  one 

interferer.  We  also  discuss  possible  generalizations  of  this  special  multiple- 

user  interference  model.  Finally,  in  Appendix  A,  we  present  a  derivation  of  a 

special  probability  function  P  (u,x)  which  plays  an  essential  role  in  this 

e 


CHAPTER  II 


HOPPING  PATTERNS  FOR  THE  FREQUENCY-HOPPED 
MULTIPLE  ACCESS  COMMUNICATION  SYSTEM 

In  a  frequency-hopped  (FH)  communication  system,  the  spreading  of  the 
spectrum  is  achieved  by  hopping  the  frequency  of  the  carrier  signal  at 
regular  intervals.  The  hopping  pattern  is  generated  by  applying  a  random 
or  pseudorandom  sequence  of  inputs  to  a  frequency  synthesizer.  Typically, 
the  available  RF  bandwidth  is  partitioned  into  q  nonoverlapping  frequency 
intervals  called  slots ,  and  the  q  different  frequencies  generated  by  the 
frequency  hopper  are  the  center  frequencies  for  these  slots. 

For  the  mobile  communication  system  described  in  [1],  some  spectral 
division  is  made  so  that  the  mobile-to-base  (upstream)  band  does  not  overlap 
the  base-to-mobile  (downstream)  band.  We  may  assume  that  the  base  is  able  to 
transmit  signals  to  the  mobiles  synchronously.  This  will  add  a  time  delay  to 
some  of  the  speech  communications,  although  the  delay  is  usually  negligible. 

We  also  assume  that  the  separation  between  two  center  frequencies,  Af  ,  is 
greater  than  1/T^,  T^  being  the  chip-time,  so  that  there  will  be  no 
interference  between  two  adjacent  frequency  slots,  and  the  FH-DPSK  system 
will  perform  like  a  frequency  division  multiplexing  system  with  the  added 
advantage  of  being  a  spread-spectrum  system. 

For  our  purposes,  a  full-hit  occurs  whenever  two  transmitters  are 
transmitting  in  the  same  frequency  slot  over  a  whole  chip-time  period,  and  a 
partial-hit  occurs  whenever  two  transmitters  are  transmitting  in  the  same 
frequency  slot  over  part  of  a  chip-time  period.  Our  objective  in  this  chapter 
is  to  produce  large  numbers  of  frequency-hopping  patterns  for  the  upstream 
and  downstream  communications  such  that  the  number  of  hits  between  two 
hopping  patterns  is  also  small. 
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2.1  A  Hopping  Pattern  Design  for  the  Asynchronous  Upstream  Communications 


In  this  case,  the  mobiles  are  free  to  transmit  signals  to  the  base 
asynchronously  so  the  number  of  hits  between  two  patterns  (with  periods  N) 
over  all  possible  time  shifts  must  be  considered.  We  may  assume  for  the 
moment  that  the  system  is  pseudoasynchr onous  (asynchronous  but  with  chip- 
synchronization)  so  that  we  can  consider  the  number  of  full-hits  between  two 
patterns.  If  x  =  (Xq,x^» . . . is  a  vector  of  length  N,  we  define  the 
left-shift  operator  T1  as  follows: 


T  (x)  (x.,xi+1,xi+2,...,xN_1,x0,...,xi_1),  0  <  i  <  N-l  (2.1 

If  2.  is  also  a  vector  of  length  N,  then  the  Hamming  distance  dH(x,^)  between 
x  and  is  defined  as 

N-l 

d  (x,£)  =  I  x(x  ,v.)  (2.2 

H  i-0  11 

where  x(j,k)  =  1  if  j  £  k  or  0  if  j  =  k  is  an  indicator  function.  It  is 
clear  that  if  d  (x,T1(y))  is  equal  to  N  for  all  i,  then  there  will  be  no 

n  —  — 

full-hits  between  x  and  v  over  all  time  shifts.  There  are  no  partial  hits 
either.  On  the  other  hand,  if  d  (x,T1(v))  =  N  and  d  (x,Ti+1(v))  =  N-l  (or 

ri  —  “  H  —  ~ 

vice  versa)  then  there  is  a  full  hit  for  one  time  delay,  no  hits  for  another 
time  delay,  and  a  partial  hit  for  any  intermediate  time  delay.  If 
dH(x,TX(£))  =  d^,(x,T^+'*'(y^) )  =  N-l  for  a  fixed  i,  then  it  is  possible  to  have 
two  partial-hits  between  x  and  y_  when  the  chip-synchronization  condition  is 
lost.  We  shall  present  a  construction  method  of  designing  a  large  class  of 
hopping  patterns  in  which  any  two  patterns  will  have  at  most  one  partial 
(or  full)-hit  when  used  asynchronously;  that  is,  if  d  (x,T1(v))  =  N-l,  then 

n  —  “ 

d^(x,T^  *■(£))  =  d^(x,T^+^ (^) )  =  N.  First,  however,  consider  the  following. 


Suppose  GF(q)  is  a  finite  field  ([7])  with  a  as  a  primitive  element. 

Solomon  ([C],[y])  had  found  a  set  of  q  Reed-Solomon  (RS)  sequences  that  can 

be  considered  as  hopping  patterns.  These  sequences  are  of  the  form 

2  \-l 

x  =  a  (1,1,. ..,1)  +  (l,a,a  ,...,a*  )  with  N  _<  q-1 ,  a  €  GF(q) ,  a  a  primitive 

X  X 

element  of  GF(q).  Thus,  the  £-th  coordinate  (entry),  0  _<  £  <_  N-l ,  of  x  is 

? 

x^  =  a^  +  a  .  The  coordinates  of  x  are  not  repeating  simply  because  if 

2,  in 

£  f  m  (0  <  £,m  <  N-l  <  q-2)  then  a  +  a  fa  +  a  thus  x.  f  x  .  Note  also 
—  —  —  x  x  £  m 

that  if  jc  belongs  to  the  set  of  sequences,  then  tSc  is  not  a  member  of  the 

1 

set.  Consider  now  two  distinct  sequences  x  and  y  with  x^  =  +  a  , 

y  =  a  +  a  ,  a  fa.  Obviously  x  f  y  for  all  £,  thus  we  have  no  hits 
&  y  x  y  jL  aj 

in  the  synchronous  case.  This  is  equivalent  to  saying  that  d  (x,v)  =  N.  But 

in  the  pseudoa synchronous  (asynchronous  system  with  chip-synchronization) 

case,  x  and  T1(y)  (i  f  0)  can  have  at  most  one  full-hit  because  the  equation 
£  i+£  i  £ 

a+a=a+a  or  (1-a  )a  =  a  -  a  has  exactly  one  solution  in  GF(q), 
x  y  x  y 

and  the  solution  corresponds  to  a  hit  only  if  £  €  {0,1 , . . . ,N-1) .  Thus  when 

N  <  q-1,  d  (x,T1(v))  =  N-l  for  some  i,  0  <  i  <  N-l;  and  when  N  =  q-1, 

H.  —  ~ 

d  (x, T1^))  =  N-l  for  all  i,  0  <  i  N-l.  Therefore,  when  N  =  q-1,  any 
two  sequences  from  the  set  will  almost  always  have  two  partial-hits  when  used 
asynchronously  (with  no  chip- synchronize tion)  .  This  is,  of  course,  highly 
undesirable. 

Example .  Consider  GF(ll)  so  that  q  =  11,  N  =  10.  Take  a  =  2 

(210  e  1  (mod  11)),  a  =  0,  a  =1.  Hence  x  =  (1,2,4,8,5,10,9,7,3,6), 

x  y  - 

Y_  =  (2,3,5,9,6,0,10,8,4,7).  It  is  easy  to  see  that  d^(x,y^))  =  N  and 
dH(x,Ti(y))  =  N-l  V  i,  1  <  i  <  9. 
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Although  in  the  example  we  chose  the  length  N  to  be  equal  to  q-1,  in 
general  the  sequences  with  length  N  <  q-1  are  more  attractive  because  Ky 
suitable  search  among  the  set  of  the  q  sequences  we  may  be  able  to  find  a 
few  with  no  hits  for  some  shifts.  However,  this  method  of  generating  RS 
sequences  suffers  a  major  disadvantage:  only  q  patterns  can  be  generated. 

In  order  to  generate  a  large  class  of  hopping  patterns  with  the  desired 
property  (maximum  one  hit  between  any  two  patterns  asynchronously),  we  may 
consider  sequences  of  length  N  <  q-1.  In  particular,  we  assume  N  is  a 
divisor  of  q-1.  To  this  end,  suppose  a  sequence  of  length  N  is  to  be 
generated : 

x  =  a  (1,1,1, . . .  ,1)  +  b(l,a  ,a  , . . . ,a  )  (2.3) 

where  W  =  (q-l)/N,  €  GF(q)  and  b€GF(q).  For  0  £  j  W-l,  define  a  collection 

of  sequences  G(cr*)  as 

G(aj)  =  (x  :  a  C  GF(q),  b  =  aj }  (2.4) 

X 

then  the  set  D  of  sequences, 

W-l 

D  =  {x  :  x  €  G(aJ)  ,  0  <_  j  <  W-l}  =  GG^)  (2.5) 

j-0 

is  a  collection  of  nonconstant  sequences  since  b  #  0.  From  our  previous 

discussion,  it  is  obvious  that  each  group,  say  G(ctJ),  contains  q  sequences 

of  nonrepeating  coordinates.  Moreover,  two  sequences  in  the  same  group 

will  have  no  full-hit  when  used  synchronously,  have  at  most  one  full-hit 

when  used  pseudoasynchronously ,  and  are  not  cyclic  shifts  of  each  other. 

Furthermore,  all  the  sequences  in  D  are  distinct  and  are  not  cyclic  shifts 

i  k 

of  each  other.  To  see  this,  let  x  €G(aJ)  and  v  €  G(a  ),  0  £  j  ,  k  <_  W-l.  Suppose 
x  =  for  some  i,  it  means  that  for  fixed  a^  and  j,  we  can  find  fixed 
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i  .  ,  j  £W  ,  k  (i+£)W  iW  j,  k+iW-j  ,  , 

a  ,k,  and  1  such  that  a  +  a  a  =  a  +  a  a  ora-a=aaJ[a  -1] 

y  x  y  x  y 

1^4-  i  W”  i 

holds  for  all  £  (a  /  a  ,  k  ^  j  so  a  ^  1) .  But  this  is  impossible, 

x  y 

Actually  since  we  take  b  (Eq.  (2.3))  to  be  nonzero  and  b  =  a*'  for  0  £  j  £  W-l 

we  automatically  exclude  all  constant  and  cyclic  sequences  from  D  because 

G(a~*  ^)  =  T^{G(a^)}  for  0  _<  j  _<  W-l,  0  _<  i  _<  N-l,  and  W  =  (q-l)/N.  Note  that 

£W  i  k+iW-i 

there  can  be  at  most  one  solution  of  £  to  the  equation  a  -a  =a  aJ[a  J-l] 

xv 

for  fixed  a^,a^  (ax  ^  a  ),  ^>1  (k  /  j)»  and  i,  so  we  can  conclude  that  the 

maximum  number  of  full-hits  between  any  two  sequences  in  D  is  1  (if  k  =  j , 

i 

x  and  v  are  in  the  same  group  so  we  are  done:  if  a  =  a  ,  since  a  -1^0 

—  —  x  y 

so  £  does  not  exist).  Hence  d^(x,T1(y))  _>  N-l  for  x,y^  G  D,  0  _<  i  _<  N-l. 

i  k 

Suppose  €  D.  Suppose  also  that  x  €  G(cr),  ^  G  G(a  )  and 
dH(x,T±(y_))  =  N-l  for  some  fixed  i,  0  _<  i  _<  N-l.  This  means  that  exactly  one 
coordinate,  say  the  £-th  coordinate  of  x  and  T^(y)  is  the  same.  Hence  if  we 
let 

1  .  n±  W  2W  (N-l)W. 

x  -  ax(l >1 »1 > • • • jl)  +  a  (l,a  ,a  ) 


/ITT  nx  k,.  W  2W  (N-l)W. 

y  ct  ( 1  $ ct  9 ot  9 •  •  •  9 ot  ) 


(2.6) 


then  we  must  have 


j  £W  ^  k  (i+£)W 

a  +  a  a  =  a  +  a  a 
x  v 


(2.7) 


i  i  w 

Now  if  it  €  G(cr  )  ,  _r  =  a(l,l,.<>,l)  +  er  (1  ,ct  ,...),  then  we  claim  that  there 
k 

is  s  €  G(a  )  such  that  d  (r,T  s)  is  also  equal  to  N-l.  In  particular, 

—  rl  —  — 

s  can  be  chosen  with  a  =  a  + (a  -a  )  such  that  r  and  T1s  will  collide  at 
—  s  r  y  x  —  — 

exactly  the  same  £-th  coordinate.  To  see  this,  let 


s  =  a  (1,1,. ..,1)  +  a  (l,ct  ,a  ,...) 
—  s 


Then  we  must  have 


j  W  ^  k  (i+Jl)W 

a  +  a  a  =  a  +  a  a 
r  s 


(2.8) 


Combining  Eqs.  (2.7)  and  (2.8),  we  can  solve  for  a^.  Thus,  we  have  proved 
the  following  proposition. 

"i  Ic  i 

Proposition  2.1.  If  x  €  G(or),  y  €  G(a  )  and  d^(x,T  (^) )  =  N-l  for  some  i, 

i  k 

0  <  i  <  N-l,  then  for  every  r_  €  G(aJ),  there  is  £  in  G(a  ) 

such  that  d  (r,T1(s))  is  also  N-l. 

H  —  — 

In  a  second  case,  suppose  x,£  are  defined  as  in  £q .  (2.6)  with  0  £  j, 

k  £  W-2 .  Then  if  £  €  G(<r'+1),  we  claim  that  there  is  £  in  G(ak+1)  such  that 

d  (r,T1(s))  is  equal  to  N-l.  Similarly,  £  can  also  be  chosen  with 

a  =  a  +  a(a  -a  )  so  that  r  and  T1(s)  will  also  collide  at  exactly  the  same 
s  it  y  x  —  —  J 

7-th  coordinate.  This  is  obvious  once  we  note  that 


_  i+1  iw  ^  k+l  (i+OW 

a  +  a  a  =  a  +  a  at 
r  s 


(2.9) 


and  solve  for  ag  together  with  Eq.  (2.7).  Observe  that  G(a^)  =  T^-’G(a  ) }.  Hence  if 
k  =  W-l,  everything  that  we  have  just  stated  applies  except  £  will  be  in  G(a^) 
and  d  (r.T'*'  ^s)  is  N-l.  If  j  =  W-l,  r  will  be  in  G(a^)  and  d  (T  ^r,T^(s))  =  N-l . 

tl  —  ri  —  — 

Corollary  2.1.  If  x  €  G(a^),  y  €  G(a^)  and  d„(x,T^(v))  =  N-l  for  some  i, 

H  —  — 

0  <  i  <  N-l,  then  for  every  r_  g  G(a^+^)  there  is  s  in 
k+l  i 

G(a  )  such  that  d^(_r,T  (s))  =  N-l.  Note  that  0  <  j, 
k  <  W-l,  and  G(aW)  =  T^Gfa0)}. 


Suppose  now  we  want  to  know  when  d  (x,T  (y))  and  du(x,T  (y))  will 

both  equal  N-l.  Intuitively,  for  a  given  field  GF(q),  this  condition  should 

depend  on  the  pattern  length  N  and  hence  on  W.  As  usual,  we  assume  that 

and  T1(y)  collide  exactly  once  at  the  £-th  coordinate.  For  d  (x,T*+^(y))  to 
i+1 

equal  N-l,  and  T  (y^)  must  then  collide  exactly  once  at  some  other 
coordinate,  say  the  m-th  coordinate,  where  m  is  necessarily  different  from  l . 
Thus  we  have  the  simultaneous  equations: 

,  j  iV  .  k  (i+£)W 

a  +aa  =  a  +  a  a 
x  y 

,  j  mW  ,  k  (i+l+m)W  (2.1 

a  +  a  a  =  a  +  a  a 

y  y 

which  reduce  to 

mW  j,  (i+l)W+k-j  ..  £W  j  .  iW+k-j  . 

a  aJ  (a  J-l)  =  a  aJ (a  J-l)  (2.1 

T  „  kl  iW+k-j  .  k2  (i+l)W+k-j  ,  ,  .  00  j  ^ 

Let  a  =  a  -1,  a  =a  -1,  and  let  a  denote  0.  Then  we 


may  solve  for  m: 


m  =  l  + 


krk2 


Since  a4  =1  and  m  must  be  an  integer,  Eq.  (2.12)  makes  sense  only  when 

(k  -k  )  mod(q-l)  is  a  multiple  of  W.  Equivalently,  w|(k  -k  )  mod(q-l). 

kL  2k? 

Note  that  in  Eq.  (2.11)  a  =0  only  when  k  =  j  and  i=0;a  =0only 

when  k  =  j  and  i  =  N-l.  But  then  x,y  will  be  in  the  same  group  and  so 


d  (x,y)  =  N  and  hence  there  is  no  need  to  consider  d  (x,T  (y))  for  i  =  0 

n  —  —  n  —  — 


and  i  =  N-l.  Thus  we  assume  that  we  do  not  consider  the  case  k. 


00 ,  or 


Proposition  2.2.  Suppose  x  €  G(er),  y  £  G(a  ).  For  each  i,  0  <  i  <  N-l, 


let  k^,k?  be  such  that 


1  iW+k-j  ,  ,  0 

a  =  a  -1  f  a 


2  (i+1 )W+k-j  ,  ,  o 

a  =  a  -1  r  a 


If  W|  (k-k  )mod(q-l)  for  some  i,  then  d  (x,T^(v))  =  du (x ,T^+^ (y) )  =  N-l 
1  Z  ri  —  —  H  —  *” 


for  these  i.  Otherwise,  d  (x,T1(y))  =  d^Cx.T1  (v))  =  N-l  can  never  happen 

H  —  —  H  —  ~ 


7  1. 


Our  problem  now  reduces  to  the  extraction  from  D  of  a  subset  D  of 

g 


sequences  which  satisfy  our  constraints.  For  a  fixed  q,  there  are  altogether 


W*q  =  (q-l)q/N  sequences  in  D.  To  pick  D  from  D,  we  do  not  have  to  generate 

g 


all  the  sequences  and  compare  them.  Indeed,  Proposition  2.2  provides  a 
reasonably  fast  and  easy  algorithm  to  find  all  the  good  sequences.  Meanwhile, 

g 

let  us  introduce  a  few  more  definitions.  We  say  a  group  G(a  )  is  a  good  group 


if  the  maximum  number  of  partial-hits  between  any  two  sequences  within  G(a  ) 


is  less  than  or  equal  to  one,  otherwise,  G(a  )  is  a  bad  group .  Likewise,  a 

r.  -  .  .  s , 


group  G(a  )  is  a  relatively  good  group  with  respect  to  G(a  )  if  the 
maximum  number  of  partial-hits  between  any  two  sequences  from  the  two  groups 

is  less  than  or  equal  to  one.  A  set  D_  of  mutually  good  groups  is  such  that 

g 

all  the  groups  inside  D  are  relatively  good  groups  with  respect  to  one 

g 

another . 

Suppose  now  we  have  a  given  field  GF(q),  a  primitive  element  a  of  the 
GF(q),  and  a  fixed  sequence  length  S,  N  divides  q-1 .  Thus  W  =  (q-l)/N  is 


determined. 


We  first  set  j=0,  k  =  0  in  Proposition  2.2.  If  W|  (k^-k^roodCq-l)  for  some 

i,  then  there  exists  x,y_€  G(a°)  such  that  d^x.T1^))  =  Dh (x ,T1+1 (y) )  =  N-l  for 

these  i.  Hence  there  will  be  two  partial-hits  between  x  and  y  when  they  are 

used  asynchronously.  Thus  G(a°)  is  a  bad  group.  Then  Corollary  2.1  implies 

that  G(a  ),G(a  ),..., G(a^  are  all  bad  groups,  hence  the  whole  set  D  is 

bad.  Therefore,  we  have  to  start  over  with  a  larger  q  or  a  smaller  N. 

Suppose  w|  (k^-k^modCq-l)  for  no  i,  then  G(a^)  is  a  good  group .  Hence  Dg 

consists  of  at  least  a  single  group  G(a^)  (or  G(a~*),  0  _<  j  £  W-l ,  but  we  shall 

always  pick  G(a^)).  Note  that  checking  the  conditions  in  Proposition  2.2  is 

equivalent  to  comparing  all  the  sequences  (since  x,y^  are  arbitrary)  and  their 

i  k. 

time  shifts  from  groups  G(aJ)  and  G(a  ),  but  in  a  much  simpler  manner. 

Next  we  set  j  =0,k=l  in  Proposition  2.2  and  check  the  stated  conditions. 
Similarly,  if  w|  (k^-k^mod^-l)  for  some  i,  then  G(aS  is  a  relatively  bad 
group  with  respect  to  G(a^).  Otherwise,  G(a'*‘)  is  a  relatively  good  group  with 
respect  to  G(a^).  Proceeding  with  j=0,  k  =  2,3 . W-l,  we  can  then  find 

all  the  relatively  bad  groups  with  respect  to  G(a^),  i.e.,the  set  of  all 

0  o  t>i  k-)  b 

relatively  bad  groups  to  G(a  )  is  RBG(a  )  =  (G(a  ),G(a  G(a  n) }  for 

some  b  ,  1  <  b  <  W-l.  For  convenience  we  shall  write  RBG(a  )  =  {b,,b„,...,b  ; 
n  —  n  —  1  2  n 

Then  Corollary  2.1  implies  that  the  set  of  relatively  bad  groups  with  respect 

to  G(a'*')  is  RBG(a^)  =  {b,+l,b_+l, . . .  ,b  +1}.  Inductively,  all  the  RBG(a'*)'s 
-  12  n 

can  be  found  without  difficulties.  Observe  that  Corollary  2.1  also  implies 
that  if  G(ar)  is  relatively  good  with  respect  to  G(aS),  then  G(ar+'*')  must 

S+l  .  t  , 

be  relatively  good  with  respect  to  G(a  ).  With  all  this  information,  it 


is  now  extremely  easy  to  construct  D  ,  the  set  of  mutually  good  groups,  by  a 

6 

simple  deletion  process. 
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Example.  Consider  GF(281)  and  let  a  =  3  be  a  primitive  element  of  GF(281). 

The  largest  possible  N  is  8,  so  W  =  35.  D  consists  of  35  x  281  =  9835 
codewords.  Each  codeword  x  is  of  the  form  x  =  a  (1,1,..., 1)  + 
b(l,  60, 228, 192, 280, 221, 53, 89) ,  where  ax  =  0 , 1 ,2 , .  .  .  ,  280 , 
b  =  1 , 3 , 9 , 27 ,81 , . .  .  ,194 ,20 ,  i.e.jb  =  . . .  a34 .  Proposition  2.2  and 

Corollary  2.1  imply 

RBG  (ct°)  =  {7,10,13,22,25,28} 

RBGCa1)  =  {8,11,14,23,26,29} 

RBG (a“)  =  {9,12,15,24,27,30} 

RBG(a3)  =  {10,13,16,25,28,31} 

RBG(a34)  =  {6,9,12,21,24,27} 

To  find  D  ,  we  first  find  the  sets 

g 

DQ  =  D  -  RBG(a°) 

=  {0,1,2,3,4,5,6,8,9,11,12,14,16,17,18,19,20,21,23,24,26, 
27,29,30,21,32,33,34} 

Dx  =  DQ  -  RBGfa1) 

=  {0,1,2,3,4,5,6,9,12,15,16,17,18,19,20,21,24,27,30,31,32,33,34} 

D,  =  D  -  RBG (a2) 


Then  D  is  just  the  last  set  D  at  the  end  of  this  process.  In  this  case, 
g  n 

the  process  ends  at  the  8-th  step,  hence  n  =  ]D  |  =9  (we  always  start  with 

g 

a  i  no/  r  /■  ori 

the  0-th  step).  Thus  =  {G(a  ),G(a  ),G(a  ),G(a  ),G(a  ),G(a  ),G(a  ),G(a"  ) 
21 

G(a  )}.  Notice  that  the  choice  of  D  is  not  unique,  but  the  number  of 

g 

sequences  in  is  fixed  (9  x  281  =  2529). 


Some  good  designs  of  D  's  with  various  q  and  N  are  exhibited  in  Table  2.1 


Notice  that  for  most  cases  N  /q  <  0.1.  It  is  conjectured  that  there  exists 

max  ^ 


a  stronger  version  of  Proposition  2.2  which  can  provide  a  deeper  understanding 


of  the  structural  properties  on  D,  hence  on  D  .  For  N/q  <  0.03,  the  codewords 

O 


in  the  D  's  of  Table  2.1  appear  to  be  excellent,  i.e.,d„(x,T  (y))  =  N-l  for  a 

g  '  n  —  — 

very  few  number  of  i’s  when  compared  with  N.  In  fact,  the  complete  Hamming 


distance  distributions  of  each  D  may  be  evaluated  through  the  help  of 

s 


Propos-  ’  :cn  2.1  and  Corollary  2.1. 


We  have  then  succeeded  in  constructing  a  design  D  of  good  hopping 

S 


patterns  in  which  any  two  patterns  will  have  at  most  one  hit  (partial  or  full) 
when  used  asynchronously. 

2.2  A  Hopping  Pattern  Design  for  the  Synchronous  Downstream  Communication 

In  this  case  we  assume  that  the  base  is  able  to  transmit  signals  to  the 
mobiles  synchronously.  This  will  add  a  time  de’ay  to  some  of  the  speech 
communications,  although  the  delay  is  usually  negligible. 

Suppose  we  have  fq  available  frequencies  in  the  downstream  band  (t  an 
integer,  q  a  prime  or  prime  power).  Partition  the  frequencies  into  t 
disjoint  groups.  For  each  group,  map  the  q  elements  of  the  finite  field 
GF (q)  isomorphically  to  the  q  distinct  frequencies.  Then  we  apply  Reed- 
Solomon  (RS)  sequences  (with  N  =  q-1)  discussed  in  Section  2.1  to  each  group. 
Hence,  for  each  group,  we  get  q  sequences  of  length  N  =  q-1,  and  there  will 
be  no  hits  between  any  two  sequences  when  used  synchronously.  Note  that  the 
method  we  mention  here  is  just  a  special  case  of  what  we  have  discussed  in 


Section  2.1. 


In  practice,  the  pattern  length  N  is  usually  a  power  of  2.  We  know  that 
17  is  a  prime,  so  patterns  of  length  N  =  16  are  obtainable.  However,  33  is 
not  a  prime,  hence  patterns  of  length  32  may  be  constructed  by  annexing  two 
patterns  of  length  16  together.  This  simple  annexing  procedure  can  certainly 
be  carried  out  for  larger  N's.  One  other  possibility  is  to  use  extended 
Reed-Solomon  sequences.  For  example,  when  q  =  33,  we  can  first  generate  32 
sequences  of  length  31,  and  then  to  the  end  of  each  of  these  sequences  we  add 
a  parity  check  element.  These  32  extended  sequences  (of  length  32)  have  no 
hits  when  used  synchronously. 


CHAPTER  III 


PROBABILITY  OF  ERROR  ANALYSIS  OF  A  NEAR-OPTIMUM  FH-DPSK  RECEIVER 
WITH  ONE  INTERFERENCE  AND  A  NONFADING  CHANNEL 

In  a  multiple-access  communication  system,  cross-user  interference 
must  be  considered  in  order  to  gain  a  true  measure  of  the  system  performance 
Since  a  bit  error  rate  performance  analysis  for  a  single-user  and  a  non¬ 
fading  channel  has  already  been  given  in  [3],  we  now  describe  a  similar 
analysis  of  the  near-optimum  FH-DPSK  receiver  with  at  most  one  interference 
at  the  n-th  time-chip  over  a  nonfading  channel. 

Suppose  that  the  system  is  chip-synchronous.  If  the  i-th  codeword  is 
being  transmitted,  the  DPSK  receiver  must  consider  two  frames  of  data  in 
order  to  make  a  decision.  Moreover,  if  the  n-th  time-chip  is  interfered 
(a  full-hit)  by  another  user  transmitting  _  the  £-th  codeword  with  the  same 
power,  then  the  receiver  input  over  the  2T-second  interval  may  be 
represented  as 

_  N-l 

r(t)  =  v/2P  1  p  (t-kT  )  [cos(y,  t+e ,)  +  5,  cos(w,t  +  ij>  )] 

kr\  -A-  C  K.  K  K.U  K.  Ill 

=0  c 
N-l 

+  /2P  I  p  (t-T-kT  )  [h.  cos(u).  t  +  6  )  +6  h  cos(m.t  +  $  )] 

,  _n  T  c  lk  k  k  kn  £m  km 


+  n(t)  (3 


where  h.,  is  the  ±1  element  of  the  i-th  row  and  k-th  column  of  the  Hadamard 
lk 

matrix,  n(t)  is  white  Gaussian  noise  (WGN)  with  two-sided  spectral  density 


NNq  that  is  thermally  generated  in  the  receiver,  and  5  the  Kronecker 

delta  function.  The  Q^'s  and  are  modelled  as  independent  identically 
distributed  (i.i.d.)  random  variables,  with  uniform  densities  over  ( 0 , 2 tt ] 
since  we  are  employing  noncoherent  phase  signaling. 

Recall  (Eq.  (1.6))  that  the  near-optimum  receiver  computes 
N-i 

A(H.)  =  Z  [<VhjkXk}  +  1  •  j  -  0,1 ... .  ,N-1 


k=0 


k  jk*k' 


and  decides  the  hypothesis  whenever  A(H^)  >A(H.)  for  j  =  0,1,...,N-1 

and  j  f  i.  Then  given  the  i-th  codeword  was  transmitted  so  that  H  is 

i 

true,  the  conditional  word  probability  that  the  receiver  makes  an  incorrect 
decision  is 


— 

N-l 

p[€ |h  ]  =  P 

U  [ A(H  )  >  A(H  ) ] 
j=0  J 

|  H . 

►  '  1 

- 

iti 

/ 

- 

Equation  (3.3)  is  bounded  from  above  by  the  union  bound,  thus 

N-l 

P[€Jh.]  <  Z  P[A(H.)  >  A(H. ) I H . ]  ( 

1  -  j=0  J  11 


Conditioning  on  H  ,  Q,  and  0  ,  and  by  using  Eq .  (3.2)  we  obtain 
Z  —  m 


P  [  A  (H . )  >  A(H .  )  I H .  ,H  ,0,4>  ] 

j  x  i  z  —  m 


=  P[  Z  {(X  -h.,X’)2+(Y  -h. ,Y’)2} 
k  ^  j  kikk  kikk 


Z  { (X  +h  X' ) 2+(Y  +h  , ,  Y ' )  2}|H.  ,H  ,0,1>  ] 
an  k  ik  k  k  lk  k  1  l  Z  —  m 


> 


( 


where  J.  *  ik:h  *  -h  /•  Here  the  terms  for  which  h  =  h  are  disregarded 
J  ik  j  k  ik  j  k 

since  they  do  not  affect  the  comparison.  Also,  J  will  contain  exactly  N/2 
elements  by  the  orthogonal  property  between  the  rows  of  the  Hadamard  matrix. 

From  [3],  we  know  that  all  the  random  variables  appeared  in  Eq.  (3.5) 
are  mutually  independent  Gaussian  r.v.'s.  By  definitions  (Eq.  (1.4)),  we 
calculate 


EI\I 

H.  , 

l ' 

'V 

,0, 

.*  ] 
m 

=  /P/2  T 

c 

[cos9,  +  5,  cosO  ] 
k  kn  m 

ElV 

H.  , 

l 

>H„ . 

A 

1  'c*'  1 

i 

m 

=  -/p/2  T 

[sin0,  +  £.  sin$  ] 
c  k  kn  m 

E[X,  1 

H. 

,0, 

i 

=  /p/2  T 

[h . ,  cos0,  +6.  h  ci 

k1 

l 

m 

c 

ik  k  kn  Jim 

EI\I 

H.  . 

l  ’ 

,0, 

,*  ] 
m 

=  -/P/2  T 

[ h . ,  sin9,  +6,  h  i 
c  ik  k  kn  Jim 

and  obtain  that  each  r.v.  has  common  conditional  variance  NT  /4 . 

o  c 

We  then  further  define  four  normalized  independent  Gaussian  r.v.'s  as 


\  ■  , 

/N.T  /2 
0  c 


.  Y.  -  h.,  y; 
„  a  k  ik  k 

*  \  ~ - 

/N  T  /  2 
0  c 


.  X,  +  h.,X' 
c  A  k  ik  k 

k  /NT  /2 
0  c 


n  AYk+hikYk 

’  Uk  / - 

/N'  T  /2 
0  c 


(3.7) 


so  that  their  respective  distributions  are 
f  n(0,l) 


n(/f^[l-h.nh^cosVl) 


l(0,l) 


n(-/PV^[l-hinhjim]sin^,l) 


n(2/PT  /N  cose,  ,1) 

C  J  (l+h.._h,_  ) 

n(2 /PT  /N„tcos6  +  - 


cos<J>_ ]  ,1) 


n(-2/PT  /N.  sinQ,  ,1) 
c  0  k 


°k~ 


_  (1+h.  h  ) 

[  /NA[sin0  +  - —  sin*  ]  ,1) 

c  0  n  2  in 


(3.8> 


where  n(p»l)  denotes  the  standard  Gaussian  r.v.  with  mean  u  and  unit  variance 
Consequently,  our  problem  in  solving  Eq.  (3.5)  is  thus  reduced  to  determining 


P  [  A  >  B|H  »H  »0,$  ]  =  /  J  f  (a,  B)dctd£  =  /  f_(3)  J  f.(a)dadB  (3.9) 
1  *  m  6=0  a=6  0  B  £  A 


2  2  2  2 

where  A  =  Z  [A,  +  B  1  and  B  =  Z  [C,  +  D,  ]  are  independent.  A 
k  €  J  “  k  k€Jk  k 

j  3 

general  solution  to  this  complex  probability  function  can  be  found  in 


Appendix  A  in  which  the  results  are  stated  in  terms  of  the  parameters 


s  =  {(EfAj2  +  (E[Bkl)2} 

and  ^ 

s’  =  l  { (E [ C  ] ) 2  +  (E[D])2} 
k  €  J  k  k 

j 


(3.10 

(3.11 


Since  the  single  interfnece  occurs  at  the  n-th  time-chip  (n  arbitrary), 

and  the  h  's  are  equally  likely  to  be  +1  or  -1,  we  have  to  consider  the 

following  three  cases  separately. 

Case  (i).  n  $  =  { k :  h . ,  =  -h.,}. 

-  j  lk  jk 

Obviously  this  case  happens  with  probability  %  since  |  J.  j  =  N/2. 

It  also  implies  that  the  interference  signal  will  not  appear  in  the 
comparison  of  Eq.  (3.5).  Hence  s  =  0,  and  s'  =  2PT/Nq.  The 
analysis  is  thus  reduced  to  the  simpler  no-interference  case  discussed 
in  [3].  In  [3]  (or  see  Appendix  A)  it  was  shown  that 


30 


P[A(H  )  > A(H  )|H  ,H  ,0,4  ]  = P [ A(H  )  > A(H . ) | H . , 0]  =p 

j  i  1a,  m  j  i  i 


-PT/2N  M/2-1  k  (PT/2NJ 

=  2  7  e  Z  2  Z  • 


k=0  £=0 


0_  /n  /  2+k- 

l !  I  k 


+k-l\ 
-£  ) 


(3.11 


Since  the  result  does  not  depend  on  0^  (or  j  or  i,  as  might  be  intuitively 
expected  because  the  rows  of  the  Hadamard  matrix  are  all  orthogonal  and 
have  equal  probabilities  of  being  transmitted),  the  upper  bound  on  the 
conditional  (on  Case  (i))  word  error  probability  are  obtained  directly: 

P[ €  I (i) ]  <  (N-l)p  (3.13 


The  relationship  between  the  conditional  (on  Case  (i))  bit  error  probability 
(=  P^Ci))  and  the  conditional  word  error  probability  is  given  in  f5]  as 


1(1)1  -srfPIS  l(i>1 


(3.14 


Substituting  Eq .  (3.13)  into  Eq.  (3.14)  we  obtain 


P,  (i)  -  p  =  P  (0  ,PT/2N  ) 
b  —  1  e  0 


(3.13; 


where  P^tu.x)  is  defined  in  Eq.  (A. 15). 


Case  (ii).  n  €  «?.,  and  h.  h  =  -1. 
- — —  j  m  £m 


Note  that  P[h.  h  =  ±  1]  =  *5.  This  case  thus  occurs  with 
m  £m 


From  Eq.  (3.10)  and  Eq .  (3.11)  we  get 


probability  h. 


IC  then  follows  easily  from  Appendix  A  that  the  conditional 
(on  Case  ( ii) )  bit  error  probability  is  given  by 


pb(11) 


) 


(3.1 


(Case  (iii) .  n  €  J.  ,  and  h.  h.  =  +1. 

-  J  in  Jim 

Similarly  this  case  also  occurs  with  probability  Although 
s  is  found  to  be  zero,  the  analysis  here  is  still  more  involved  since 


2PT 


4PT 


m  +  “jf —  [2cos0  cos$  +  cos  $  +2  sinG  sind)  +  sin^<i>  1 

nq  n  m  m  n  m  m 


2PT  4PTc 

— n—  +  -n — [1  +  2cosij>  ] 
N0  N0  nm 


(3.1 


where  "  0n  "  modulo  2tt,  is  not  a  constant.  Observe  that 

and  are  i.i.d.  r.  v.'s  with  uniform  distribution  over  (0,2tt], 

we  can  therefore  show  that  ¥  is  also  uniform  in  (0,2tt1.  With 

nm 

this  information  in  mind,  and  by  Fact  A. 2  in  Appendix  A,  we  thus 

obtain  the  conditional  (on  Case  (iii))  bit  error  probability  by 

removing  the  conditioning  on  'V  : 

nm 


(3.1 


Combining  all  the  results  from  these  three  cases,  we  thus  successfully 


provide  an  upper  bound  of  the  bit-error  probability  for  the  chip-synchronous 
system  with  single  interference  as 


pb  1  *5  pb(i)  +  *  pb(i:L)  +  *  Pb(i:Li) 


(3.20) 


where  P  (i) ,  P  (ii) ,  and  P,  (iii)  are  given  bv  Eqs .  (3.15),  (3.17),  and 
b  b  b 

(3.19),  respectively. 

Although  it  took  some  effort  to  have  arrived  at  Eq.  (3.20),  the  upper 
bound  is  still  not  general.  In  fact,  in  a  realistic  system,  chip-synchronization 
between  the  transmitters  is  almost  unrealizable.  As  a  consequence,  the 
receiver  will  experience  a  single  partial-hit  (we  avoid  the  possibility  that 
two  or  more  partial-hits  may  actually  happen  by  time  shifting,  also  see  Chap¬ 
ter  II)  rather  than  a  single  full-hit  at  the  n-th  time-chip.  This  problem 
may  be  remedied  by  annexing  an  independent  r.v.,I,  with  uniform  distribution 


over  (0,1)  to  the  interference  signal.  Thus  Eq.  (3.8)  becomes 


n(0,l) 


n (a/PT  /N  [1-h  h  Jcosi  ,1) 
c  U  in  Jem  m 


n(0,l) 

n(-a/PT  /N  [1-h  h  jsin|  ,1) 
c  0  in  m  m 


h(2/PT  /N_  cos6,  ,1) 

:  J; _ _  k  (i+h.  h  ) 

k  ]  n (  2/PT  /N.  [cose+G - --n  H  cosp  ]  ,1) 

l  c  U  n  £.  m 


h(-2/PT  /N  sin9  ,1) 

),  c  U  (1+h.  h  ) 

K  n(-2/PT  /N  [sine  +  a - simt1  ]  ,1) 

L  c  U  n  l  m 


(3.21) 


where  a  €(0,1). 


Conditioning  on  E,  the  conditional  bit-error  probabilities  for  a 
partial-hit  are  easily  found  to  be 


Pb(i|E  =  o)  <  P  (0.PT/2N  ) 


(3.22) 


Pb(iij  Z  -  o)  <  Pe(PTca  /Nq ,PT/2Nq) 


P,  (iii  I  Z  =  a)  <  h  P 
b  '  —  e 


(3.23) 


(3.24) 


which  are  similar  to  Eqs.  (3.15),  (3.17),  and  (3.19).  Note  that 
2 

?e (PTc°  /NQ,PT/2Nn)  <  P^(PT  o/Nn ,PT/2Nn)  which  is  convex  in  a. 


0  -  e  c  0 


P  (0.PT/2N  +  PT  /N  (a  +2o))  is  convex  and  monotone  decreasing  in  a,  and 

e  0  c  0 

2 

Pe(0,PT/2Ng  +  PT^/NgCcj  -2a))  is  monotone  increasing  in  a  (see  Appendix  A, 
Fact  A.  3)  V  0“  (0,1),  we  thus  obtain 


P,  (i)  <  P  (0 ,PT/2N  ) 
b  —  e  U 


(3.25) 


P.  (ii)  <  is  P  (0.PT/2N.)  +h  P  (PT  /N  , PT/2N  ) 
b  —  e  U  ecO  0 


(3.26) 


P  (iii)  <  k  P  (0.PT/2X  )  +k  P  (0 ,PT/2N  +  3PT  /N.)+'i  P  (0,PT/2Nn-PT  /N.) 
b  —  e  0  e  OcO  e  OcO 

(3.27) 

by  removing  the  conditioning  on  E.  Finally,  the  average  bit-error  probability 
for  the  chip-asynchronous  system  with  single  partial  interference  is  bounded 
by  Eq.  (3.20) 


P,  <  4  P,  (i)  +  k  P,  (ii)  +  k  P,  (iii) 


34 


For  the  chip-synchronous  system  with  single  full  interference,  the  upper 
bound  (Eq.  (3.20))  is 

pb  1  k  pb(i)  +  k  pb(ii)  +  k  pb  (iii) 

<  k  p  (0 ,PT/2Nn)  +  k  ?  (PT  /N..PT/2N.)  +  1/8  P  (0,PT/2Nn  +  3PT  /N_) 

—  e  U  ecu  U  e  U  cU 

+  1/8  P  (0 ,PT/2N  -  PT  /N  )  (3.29) 

e  0  c  0 

On  the  other  hand,  if  the  receiver  encounters  no  interference,  the  upper  bound 
[2]  is  just 

P,  <  P  (0  ,PT/2N  )  (  3.30) 

b  —  e  0 

These  bounds  were  evaluated  in  terms  of  the  bit  SNR,  E  /N  ,  for  several 

b  0 

values  of  N,  and  are  shown  graphically  in  Figures  3.1-3. 7.  Note  that 
PT/S0  -  S  PTc/S()  -  n  E^. 

As  one  might  expect,  the  upper  bound  is  largest  in  case  of  a  full-hit 
for  all  N's,  since  the  receiver  encounters  the  most  interference  energy. 

Figures  3.3  and  3.4  show  that  the  anti-interference  capability  of  the  receiver 
increases  drastically  for  larger  N's.  This  is  also  intuitively  expected 
since  more  signal  energy  is  present  relative  to  the  interference  signal 
energy,  for  larger  N's.  However,  there  is  a  trade-off  between  the  code  lengths 
and  the  receiver  performance  because  receiver  complexity  also  increases 
significantly  with  N. 

-3 

At  a  bit  error  rate  of  10  ,  Figures  3.3,  3.6,  and  3.7  indicate  that  a 

code  length  of  N  =  32  is  optimum  for  the  receiver  with  no-hit,  single  partial- 

hit,  and  single  full-hit, respec tive Ly .  Moreover,  for  N  =  32,  bit  SNR's  of 

-3 

about  8.2  dB,  8.3  dB,  and  8.4  dB  are  adequate  for  a  bit-error  rate  of  10 
for  the  three  cases,  respectively . 


Figure  3.3.  Upper  bounds  of  average  bit-error  probabilities  for  the  near 
optimum  receiver  with  a  single  interference  and  nonfading 
channel.  N  =  32. 
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Bit  SNR,  Ejj/Nq,  dB 

Figure  3.4.  Upper  bounds  of  average  bit-error  probabilities  for  the  near 
optimum  receiver  with  a  single  interference  and  nonfading 
channel.  N  *  64. 
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CHAPTER  IV 


PROBABILITY  OF  ERROR  ANALYSIS  OF  A  NEAR-OPTIMUM  FH-DPSK  RECEIVER 
WITH  MULTIPLE-USER  INTERFERENCE  AND  A  NONFADING  CHANNEL 

We  are  now  in  a  position  to  evaluate  the  performance  of  the  FH-DPSK 
receiver  in  the  presence  of  multiple-user  interference  over  a  nonfading 
channel.  In  Chapter  III  we  have  shown  that  the  effect  of  a  single  full-hit 
on  the  error  probability  is  always  more  severe  than  that  of  a  single  partial-hit 
Therefore,  for  the  sake  of  simplicity,  we  assume  that  the  system  is  chip- 
synchronous.  Under  this  assumption,  the  upper  bound  on  the  error  probability 
will  be  looser  than  the  one  for  a  realistic  asynchronous  system.  However, 
the  simpler  upper  bound,  which  is  a  union  bound  on  the  worst-case  error 
probability,  will  also  be  applicable  to  the  more  realistic  system. 

Suppose  the  i-th  codeword  is  being  transmitted.  Suppose  also  that  the 
n^  .n^ , . . .  ,nj-th  (0  _<  J  _<  N)  time-chips  have  interference,  each  from  a  single 
distinct  interferer,  and  the  J  distinct  interferers  are  transmitting  the 
£, £9 ,...,£  -th  (the  £.'s  may  not  be  distinct)  codewords  respectively  with  the 
same  power.  Then  the  receiver  input  over  the  2T-second  interval  is  (Eq.  (3.1)) 

_  N-l  J 

r(t)  =  /ZP  I  p  (t-kT  )[cos(u,  t+0  )  +  Z  6  cos(to  t+<j>  )] 

<  «  x  c  k  k  ^  ten  tc  m 


+  ^  APT  <t-I_kIc)  [hikcos(“kt+ek>  +  Via,  m  c°s(V+V  )! 

k=0  C  T  =  1  TXT  T 


+  n(t) 


(4.1) 


where  h  is  the  ±  element  of  the  £  -th  row  and  the 
£  m  x 

x  x 

Hadamard  matrix,  ©  ~  U(0,2tt]  for  all  k,  $  ~  U(0,2tt] 


m  -th  column  of  the 
x 

for  all  m  . 

T 


Thus,  given  H^,  3.,  (n^  ,  H  ,  Che  conditional  expectations  (Eq.  (4.1)) 

are 


EI\IVe, V 

E[Yk|H1,a,(nT}^1>Hi,^] 


/  /P/2  T  cos 9, 
c  k 


j  k  ^  9n^  9  •  •  •  »r 


y?/2  t  t  cos9  +cos<j>  ]  ,  k  =  n 


c  n^ 


/P/2  T  [cos9  +cos<)>  ]  ,  k  =  n. 


C  nJ  mJ 


/p/2  T  [cosQ,  +  E  5,  cos4)  ] 

c  k  .  kn  m 
T=1  x  x 


-/P/2  T  sm0.+  E  <S  sin<{>  ] 
c  k  ,  kn  m 
r=l  x  x 


/P/2  T  [h.,cos0  +  E  6,  h  cos<ji  ] 
c  ik  k  T=1  knT  iTmT  mT 


-/F/T  T  [ h . ,  s in 0,  +  E  6,  h  sin<J>  ]  (4.2) 

c  ik  k  ,  kn  £  m  m 

X— 1  XXX  X 


and  the  conditional  variances  are  all  equal  to  a  =  N_T  /4 .  After  normaliza- 

0  c 


tion,  it  follows  from  definitions  (Eq.  (3.7))  that 


X  -h.,  X/ 
k  ik  k 


\  = 


Wk 
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Vhikxk 
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~  n 


The  set  J.  is  defined  to  be  {k:h..  =  -h  }.  Using  Eqs.  (3.9),  (3.12),  and 
J  lk  jk 

(3.15)  we  see  that  the  conditional  bit  error  probability  is  bounded  by 


^  •  P(A>  B  H  ,0,ln  }J  H  ,*  )  ?(  Z  (A  ) ~+(B  )  >  1  (C  ) “+(D,  )“  H .....  , 

2  1  -  X  x  -  -m  2  I  k  k  k(ZJ  k  k  1 

'  j  j 

(4. 

If  n^  $  J  ,  the  interference  signal  at  the  n^-th  time-chip  will  not  affect 
the  overall  comparison.  Suppose  now  that  £  of  the  J  n  's  lie  inside  J.  . 
Since  jw  ]  =  N/2,  and  because  we  have  restricted  the  number  of  interferers 
at  each  time-chip  period  to  be  at  most  one,  it  is  clear  that 
(J-N/2)  u(J-N/2)  _<  e  _<  min(J,N/2)  where  u(n)  =  1  if  n  >  0,  0  if  n  <  0  is  the 
unit  step  function.  Thus,  the  probability  that  c  of  the  J  n  's  lie  inside  J) . 


P(E  =  e)  = 


I  ^ 

6  V 


In  this  case,  Eq.  (4.4)  suggests  that 

2 

A  =  2(N/2  -  e)  random  variables  (r.v.'s)  ~  n  (0,1) 


+  e  r.v. 's 


( 1-h .  h„  )cos<j> 
N„  in  i  m  it 

0  TIT 


+  e  r.v.'s  ~n  f-~\  — rr~  (1-h.  h„  )sin<j>  ,11,  for  some  e  of  x 's 

V  V  ho  “t  lA  J 


B  =  (N/2  -  e)  r.v.'s  ~  ^  cos\’^ 


+  (N/2  -  t)  r.v.'s 


—  2  sine  ,1 

No 


+  E  r.v.'s  ~n2r\/— rr-(2  cose  +(l+h._  h 


„  _  ) c os d  ),l) 


+  e  r.v. 's  ~  n 


— -  (2  sin6  +(.-(-h.  h  )sin<j)  ),1 
N_  n  in  l  m  m 

0  t  t  T  T 


for  some  (N/2--:)  of  k's  and  some  e  of  x's 


(4.6) 


where  n  (y,l)  denotes  the  square  of  a  Gaussian  r.v.  which  has  mean  y  and 

variance  1.  Rearranging  the  indexes  on  h.  ,  h„  ,  and  6  ,  we  have  from 

in  £  m  m 

TXT  X 

Eqs.  (3.10)  and  (3.11)  the  parameters 


PT  e 

s  =  2  — -£■  E  (1-h.  h  ) 
N-  ,  m  l  m 
0  x=l  x  x  r 


(4.7) 


PT  PT  e  PT  e 

s'  =  4  -sr  *  t  +2~ir  E  (1+h-  )  +  4  i  d+h.  h  )-coS(e  -<p  ) 

NA  2  N_  .  in  £  m  NA  ,  in  £  m  n  m 

0  0  x=l  xxx  0  x=l  xxx  xx 

Note  that  (1-h.  h  )  can  either  be  0  or  2.  Thus,  assuming  £  of  the  e 
mum 

X  XX 

(1-h.  h„  )  terms  in  Eq.  (4.7)  are  equal  to  zero,  we  obtain 
m  £  m  n 

xxx 


s  =  4  ~  (e  -  c) 
0 


PT  PT  c 

s'  =  2  +  4  IT  ;  +  8  z  COS^x 

0  0  0  x=l 


(4.8) 


where  V  =  0  -  0  (mod  2ir)  is  a  r.v.  ~  U(0,2tt].  The  solution  to  Eq.  (4.4) 

x  n  m 
x  x 

is  then  given  by  Eq.  (A. 15)  in  Appendix  A  to  be  Pe(u,x)  with  u  =  s/4  and 
x  =  s'/4. 

Fortunately,  all  the  random  variables  in  Eq.  (4.4)  are  mutually  independent. 
This  fact  enables  us  to  remove  the  conditionings  one  by  one  via  taking 
independent  expectations  with  respect  to  each  individual  random  variable. 

Using  Fact  A. 2  in  Appendix  A  successively  z  times,  we  can  remove  the 


'  V..  ,  .  .V— .a.V-.  v  v  V  .'w  V_.  ^  i  .*%— .  A.  1a. 
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conditionings  on  4^  s  and  obtain  an  upper  bound  of  the  bit-error  probability 


to  be 


^  l  r  r  f  PTc  PT  PPc 

2  (y)c  C)PQ  (  -Trie  -  c),  +  ~(3C 

o=0  2  a  6  \  N)  2N0  N0 


-  4a)j 


(4.9) 


But  P (Z  =  £)  is  binomially  distributed,  hence,  removing  the  conditioning  on  Z, 
the  upper  bound  is 


£  £  /  PT  PT 

-  p  4)5(')pe(  -#(E  -  ?>.?T  +  -  4c 

c=0  5  “  a=0  a  V  0  0  0 


(4.10) 


Finally,  removing  the  conditioning  on  E  (the  number  of  n^'s  that  lie  inside 


J.)  yields 


e  £  /  PT  PT  \ 

P  <  r  P(E-c)  E  <VV  I  <7)C(J)Pe  (  -fke  -  5>.§-  +  ir(3£  -  4o))(4.U) 
D  e  £=0  5  2  a=0  “  a  ®  \  1 0  0  0  J 


where  P(E=e)  is  given  by  Eq.  (4.5).  Thus,  given  that  J  interferers  are  present 
and  each  of  the  interferers  is  interfering  a  distinct  time-chip  period,  the 
bit-error  probability  is  bounded  by  Eq,  (4.11).  It  is  not  difficult  to  see 
that  Eq.  (4.11)  agrees  with  Eq.  (3.30)  and  Eq.  (3.29)  for  the  J  =  0  and  J  =  1 
cases,  respectively  . 

The  expression  in  Eq.  (4.11)  can  be  determined  using  a  digital  computer 


as  a  function  of  the  bit  SNR,  E,  /N  for  various  values  of  N  and  J.  Note  that 

b  U 


PT 

N  c  PT 


VN,= - —  =  — — — ,  and  N  =  2  .  Several  upper  bounds  are  displayed  graphically 

0  n  Nq  nNg 

in  Figures  4.1  to  4,4.  For  N  =  8  and  N  =  16,  the  receiver  apparently  fails 


to  operate  within  practical  limits  of  bit  SNRs  if  each  of  the  time-chips  is 
interfered  by  a  different  interferer.  For  larger  N's,  the  receiver's  anti¬ 
interference  capability  is  seen  to  be  quite  satisfactory.  In  fact,  at 


Upper  bounds  of  average  bit-error  probabilities  for  the  near 
optimum  receiver  with  single  interference  at  J  distinct  time 
chips  and  nonfading  channel.  N  =  8. 
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=  10  and  assuming  that  1/8  of  the  N  time-chips  are  each  interfered  by  a 
distinct  interferer,  the  degradation  in  bit  SNR  is  just  1.25  dB  for  N  =  32 
and  0.8  dB  for  N  =  64.  By  using  the  frequency-hopping  patterns  that  we 
developed  in  Chapter  II,  it  is  believed  that  the  number  J  of  interferers 
(to  the  receiver)  will  be  small  if  the  number  M  of  simultaneous  users  in  the 
system  is  not  too  large.  Although  multiple  interference  at  the  same  time- 
chip  period  are  possible  in  a  realistic  system,  it  will  happen  with  a  small 
probability  since  any  two  sequences  have  at  most  one  hit  asynchronously  and 
each  sequence  is  hopping  on  a  small  portion  of  the  available  frequency 
spectrum  (for  example  N/q«»  0.03  in  Chapter  II). 

For  a  fixed  number  M  of  simultaneous  users,  we  may  assume  that  each  group 
in  a  good  design  D  (Table  3.1)  accommodates  approximately  an  equal  number  of 

O 

users.  With  the  help  of  Corollary  2.1,  it  should  not  be  difficult  to  compute 
the  entire  hitting  characteristic  of  the  system  with  M  users.  Then  we  can 
deduce  the  number  of  interferers  (to  the  receiver)  and  their  distributions 

so  that  the  bit-error  probability  rates  can  be  calculated.  Typically,  N  is 

-3 

32  and  P  is  10  .  Thus,  within  practical  limits  of  the  bit  SNRs  (10  ~  20  dB) , 

the  maximum  number  of  simultaneous  users  in  the  system  can  be  found. 

As  a  preliminary  investigation,  we  may  assume  that  about  1/10  of  the 

total  |D  ]  potential  users  are  actually  transmitting  signals.  Note  that  q 
§ 

distinct  frequencies  are  available  (see  Chapter  II)  and  the  sequence  length 

is  N.  Hence,  any  frequency  f  will  appear  in  some  |D  | -N/lOq  distinct 

sequences  and  the  maximum  number  of  interferers  can  be  as  large  as 

(Id  I  *N/10q-l)N.  But  among  the  ID  1/10  users’  sequences,  some  of  them  mav 
S'  g 
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never  have  hits  with  the  receiver's  sequence  because  two  sequences' 


frequency  slots  can  be  disjoint.  On  the  other  hand,  some  of  the  sequences 


may  have  more  than  one  frequency  in  common  with  the  receiver's  sequence,  so 


that  the  number  of  interferers  can  be  much  smaller.  Furthermore,  there  are  N 


possible  time  shifts  of  a  sequence  so  it  is  very  unlikely  that  all  the 


potential  interferers  will  actually  interfere  with  the  received  signals.  It 


is  then  reasonable  to  assume  that  the  number  of  interferers  J  is  about 


ID  I •N/10q-l.  For  example,  from  Table  2.1,  we  mav  take  the  design  D 
1  g‘  '  S 


GF(577)  with  N  =  32.  |D^|  is  given  to  be  2  x  577  =  1154.  Assuming  that 


M  115  simultaneous  users  are  present  in  the  system,  then  J  is  about  5 . 


We  shall  allow  multiple  interferences  at  some  time-chip  periods.  Consequently, 


the  ways  that  J  interferers  will  distribute  in  N  time-chips  will  have  a 


Bose-Einstein  distribution.  The  bit-error  rates  can  then  be  estimated  easily 


by  using  the  function  P^Cu.x).  Although  this  kind  of  interference  analysis 


will  inevitably  expend  a  substantial  amount  of  computer  CPU  time,  it  is 


certainly  a  necessary  follow-up  to  this  thesis. 


-  •-  J 
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CHAPTER  V 
CONCLUSION 

In  this  thesis  we  have  discussed  a  design  of  short  Reed-Solomon  sequences 
in  which  any  two  sequences  will  have  at  most  one  full-l^it  or  one  partial-hit 
when  used  asynchronously.  The  result  showed  that  these  RS  sequences  have 
excellent  correlation  properties  and  therefore  are  very  suitable  for 
frequency-hopping  application,  since  they  can  extend  the  number  of  multiple- 
access  users.  We  also  discussed  the  behavior  of  the  FH-DPCK  near-optimum 
receiver  with  multiple-user  interference  in  a  nonfading  channel,  assuming 
that  each  time-chip  has  at  most  one  interferer  of  the  same  power.  The 
receiver  performance  was  found  to  degrade  gradually  with  an  increase  in  the 
number  of  interfering  signals.  Furthermore,  a  longer  codeword  length  was 
shown  to  be  more  preferable.  Finally,  generalization  of  the  multiple-user 
interference  model  was  considered. 


APPENDIX  A 


THE  ERROR  PROBABILITY  P  (u,x) 

e 

Our  point  of  departure  is  the  probability  expression 


CO  OO  OO  CO 

P[A  >  B]  =  /  /  f,  (a, 8)dadB  =  /  fR(8)/  f,(a)dadB  (A.l) 
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where  A  is  independent  of  B.  Now  suppose  that  A  and  B  are  both  being  a  sum 
of  N  independent,  squared  Gaussian  random  variables  with  nonzero  means  and  unit 
variances.  Then  A  and  B  will  have  the  noncentral  chi-squared  density 
with  N  degrees  of  freedom  given  in  [10]  as 


and 


N-2  -  (a+s) 

fA(c0  “  I  (f}  4  e  2  a>° 

2 

N-2  - ( 8+s 1 ) 

£B(B)-i(A)4  e  2  IN2(v^-).  B  >  0 
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(A. 2) 


(A.  3) 


where  s  and  s'  are  sums  of  the  squares  of  the  means  of  the  Gaussian  r.v.'s  in 
A  and  B, respectively .  The  inner  integral  in  Eq .  (A.l)  is  then  found  [10]  to  be 


/  fA(q)da  = 


where 
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is  a  generalization  of  the  Marcum  Q-function.  Observing  that  Q(/s", /ff)  is 
given  by  [10]  as 
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where  I^CvTb)  is  the  r-th  order  modified  Bessel  function  of  the  first  kind 
given  by 


m) 
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we  obtain 
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Upon  substitution  for  Eq.  (A. 8)  and  Eq.  (A. 3)  into  Eq.  (A.l)  and  by  using 
the  fact  that  f  ( 3 )  is  a  probability  density  function  we  get 
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where  v  =  s'/2,  and  y  =  2B.  The  integral  in  the  right-hand  side  of  this 


expression  may  be  solved  with  the  least  amount  of  difficulty  through  indirect 
computation.  We  first  notice  that  the  characteristic  function  of  a  positive 
random  variable  Y  may  be  put  into  the  following  form  by  using  the  power  series 
expansion  of  the  exponential  function: 


OO  00  k 

*Y(ju)  -  E{ej“Y}  *  E  (ju))k[J  |r  fY(y)dy] 

k=0  0 
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The  integral  in  brackets  is  observed  to  be  the  integral  of  interest,  where 


f^(y)  is  the  noncentral  chi-square  density  encountered  earlier  in  Eq.  (A. 2). 

Thus,  the  value  of  the  integral  is  simply  the  value  of  the  coefficient  of 

(jw)  in  the  series  expansion  of  the  characteristic  function  corresponding 

to  the  noncentral  chi-square  density  with  N  degrees  of  freedom.  That 

characteristic  function  is  given  in  [10]  as  ^(ju)  =  (l-2jw)  exp[vjco/(l-2jw)  ] 

This  function  may  be  expanded  by  first  using  the  power  series  expansion  for  the 

exponential  function,  indexed  by  .  Then  hv  using  the  binomial  expansion 
-(N/2+Jc,) 

for  each  (l-2jo>)  term,  indexed  by  m,  we  obtain 


)  =  I  I  ( N/2+JWm-A  (2ju) 


£=0 


U 


k=0  \ 


m 


\  (2j. 


kra 


(A. 11) 


^  ^  | 
where  (,  )  denotes  the  usual  binomial  coefficient  N •, 

b  b  !  ( a-b )  ! 


Finally ,  by 


.  k+r 


collecting  all  the  contributions  to  the  coefficient  of  (jw)  ,  the  value  of 
the  integral  is  found  to  be 


,k+r  ktr  (v/2)£  ( N/2+k+r 


£= 0 


V. 


V 


k+r 


+r-l\ 

-i  I 


(A. 12) 


Thus , 


N 


P[A  >  B]  =  1-2  e 


y_ 

2 


\T/0  ,  „2r+3k  , 

r=N/2  k=0  2  k! 


2k+r  k+r(v/2)* 


Z=0 


% ! 


/  N/2+k+r-l' 
l  k+r- i 


=  1-2 


N 
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2  -2u  -x 
e  e  l 

y-=M  /? 


,-r  k 
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ktr  A  /  N/2+k+r- 

0_0  ll  l  k+r 


-0 


(A. 13) 


where  u  =  s/4  and  x  =  v/2  =  s  /4.  Hence  Eq .  (A.l)  is  solved.  Notice  that  if 
s  =  0,  then  Eq.  (A. 13)  simplifies  to 


2-N/2e-x 


N/2-1 
E  2 
r=0 


-r 


r 

E 

1=0 


xV  N/2+r-A 

*  •  \  r — J 


(A. 14 


which  is  exactly  the  special  case  error  probability  expression  p  given  in  [J]. 
We  then  define 


P  (u,x)  =  ^  P[A > B]  (A. 15 

e  2 

where  P[A>B]  is  given  in  (A. 13).  It  is  obvious  that  P^Cu.x)  is  continuous 

and  differentiable  in  u  and  x.  Furthermore,  it  can  be  shown  that  P  (u,x)  is 

e 

convex  [11]  in  both  u  and  x,  monotone  increasing  in  u,  and  monotone  decreasing 

2  2 
o  a  9 

in  x  by  observing  that  — -  P  (u,x)  21  0.  — J  pe^u»x^  1.  °»  ^  Pe(u,x)  >  0, 

3u^  6  3x  6  U  e 

'b 

and  —  P  (u,x)  <  0  V  u,  x  >  0,  respectively.  Figure  A.l  and  Figure  A. 2  show 
dX  e 

graphs  of  P  (u,x)  for  interesting  values  of  u  and  x.  In  particular, 

®  E,  PT  EL  \ 

P  /u,x  =  —A  and  P  [  u  =  — rr— •  =  ^  TJ—  »  x  ]  are  exhibited  in  the  figure 

eV  2No  2  N0 J  e\  No  N  N0  / 

with  N  =  2n  =  32. 

Finally,  we  present  some  useful  facts  which  follow  readily  from  the 

convexitv  of  P  (u,x). 

e 

Fact  A.l  Given  a  >  0,  b  >  0,  and  E  ~  U(0,1), 

E{P  (u,X  =  a  +  bE)}  <  ^{P  (u,a+b)  +  P  (u,a)}  (A. 16 

e  —  2  e  e 

Fact  A.  2  Giver,  o  >  0,  a-b  >  0,  and  U(0,2tt  ], 


E{P  (u,X  =  a  +  b  cos$)  }  <  -jkp  (u,a+b)  +  P  (u,a-b)} 


(A. 17 


Fact  A.  3  Given  a  €  (0,1),  P  (0  +  — —  (a  +  2a))  is  convex  and  monotone 

e  2Nq  Nq 

PT  PTc  2 

decreasing  in  a,  and  P  (0,  -r— —  +  — —  (a  -  2a))  is  monotone 

e  2Nq  Nq 

increasing  in  a. 

Fact  A. 2  is  trivial  once  we  observe  that  g(i|j)  =  Pe(u,a+b  cosip)  is  convex  and 
increasing  in  ( 0 , tt ]  ,  and  convex  and  decreasing  in  (it,2tt]  .  Since  $  is  uniform, 
E{g(^>)}  will  simply  be  equal  to  1  / 2 tt  times  the  area  under  g(ip)  .  The  result 
then  follows  from  convexity  after  approximating  g(ip)  by  two  straight  lines. 

Fact  A. 3  may  be  verified  by  Figure  A. 3  and  Figure  A. 4. 

This  concludes  our  discussion  on  P  (u,x) . 


Figure  A. 4.  Weighted  graphs  of  c  •  P  (0,-^— —  +  — — 

e  2Nq  Nq 


(ct  -2<j))  with 


EL  PT  E 

PT/N  =  a  ^  a  =  5,  N  =  32. 

0  0  0 


Graph  1  E^/NQ  -  2  dB,  c  -  10 

2  E^/Nq  =  4  dB,  c  =  1 

3  E^/Ng  =  6  dB,  c  *  1 

4  E^/Nq  =■  8  dB ,  c  *  IQ2 

5  E^/Nq  -  10  dB,  c  =  104 

6  E^/Nq  =  12  dB,  c  =  109 
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